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Abstract. The purpose of this article is to provide new error estimates 
for a popular type of SBF approximation on the sphere: approximating 
by linear combinations of Green's functions of polyharmonic differential 
operators. We show that the Lp approximation order for this kind of 
approximation is o for functions having Lp smoothness a (for o up to 
the order of the underlying differential operator, just as in univariate 
spline theory). This improves previous error estimates, which penalized 
the approximation order when measuring error in Lp, p > 2 and held 
only in a restrictive setting when measuring error in Lp, p < 2. 

1. Introduction 

Spherical basis functions (or SBFs) have been used with much success 
in multivariate approximation theory, statistics and a multitude of other 
scientific disciplines. At the heart of the SBF methodology is the creation 
of an approximant 

by taking a linear combination of rotations of a fixed kernel {x, a) ^ k{x ■ a) 
(known as an SBF, or, sometimes, a zonal kernel). 

The success of the SBF methodology derives from its ability to generate 
approximants from data having arbitrary geometry - a desirable quality on 
spheres, where geometry of data is always essentially unstructured: for an 
arbitrary spacing, there are no regular distributions of points on the sphere, 
meaning that approximation techniques requiring grids, regular triangula- 
tions, or other geometrical props do not work in this setting. Interpolation 
[H, [3], [9] and [7] and other SBF approximation methods [II], [H], [l8], [TO] 
(see bibliography in [3j for even more examples), are both frequently used 
to fit scattered data on the sphere. 

Our focus is not on how to treat spherical data, but how to approximate 
smooth functions using SBF approximants, having access to as much infor- 
mation about the target function as necessary. The choice of coefficients 
(^^)ggH is a crucial element in the performance of the approximation, but, 
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at the outset, we are not focused on a specific method of choosing coef- 
ficients. Instead, we wish to investigate the approximation power of this 
methodology for a robust family - the polyharmonic kernels (see Definition 
13. 1[) - rather than any specific implementation or algorithm; the main con- 
cern is to establish accurate error analysis for approximation from spaces of 
polyharmonic SBFs, S{k, H) := span^g^ ^) -1-11, where a low dimensional 
space of elementary functions, 11, may be added to the span of the SBF. 

The method for gauging the approximation power is the Lp approxima- 
tion order, which measures the decay of the error in approximating from 
S{k, H) as H becomes dense in S*^. For target functions / from a class 9", 
the approximation order is the largest exponent s so that 

II/-s/,h||l,(s<^) = 0(/^^) 

where h, the 'fill distance', measures the density of H in S*^ (see the following 
section for a precise definition of fill distance). In this setting, the rate is 
given in terms of the density of the centers H, and depends strongly on the 
class 3" of target functions. 

When approximants are chosen from a predetermined linear space, inde- 
pendent of the target function, as is the case here (in contrast to nonlinear 
approximation, where the set of centers H could be chosen independent of 
/), precise approximation theory ties the Lp approximation order to the 
Lp smoothness of the target function, e.g., by measuring the error in terms 
of an Lp modulus of smoothness or by selecting target functions in an Lp 
Sobolev or Besov space, 9" = Wp or Bp^^. 

The prevailing method for estimating error for SBF (and, more generally, 
kernel) approximation has been to assume the target function resides in a 
reproducing kernel Hilbert space, often called the native space, for which 
the SBF acts as the reproducing kernel. Quite often, the native space is ac- 
tually an L2 Sobolev space. One drawback of this approach has been that it 
precludes finding faster rates for functions with more smoothness, or slower 
rates for less smoothness. Another drawback is that the Lp approximation 
orders degrade as p aberrates from 2 - see, e.g., [6] for an example of this 
criticism for 'radial basis functions' (or RBFs) in domains in M^. We remark 
that [ini Corollary 3.5], and [9l Corollary 3 (a)] are examples of this phenom- 
enon, but we place special emphasis on the results of Hubbert and Morton 
[H Theorem 3.4, 3.8], because their results are the current state-of-the-art 
for the setting of this article. We paraphrase their result. 

Theorem (Hubbert, Morton). For an SBF k having native space W^{S'^), 
and for sufficiently dense centers H, if f ^ W^{E>'^) then the SBF interpolant 
Sf satisfies: 

If f e ^|'"(S'^), Sf satisfies: 
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When p > 2, each approximation order is penahzed by subtracting a 
positive term: f — ^; when p < 2, the space of target functions is an L2 
Sobolev space of the form W^(E>'^) which is embedded in W^{S'^). This 
should be contrasted with M^^ order univariate sphne approximation, which 
provides Lp approximation order a for functions in Ty^(§'^) for a range of 
< cr < M), where M is the 'saturation' order - the rate beyond which any 
increase in smoothness fails to produce an increased rate of convergence. Our 
main results show (in Theorem 16.11 and its corollaries) for a polyharmonic 
kernel, k, satisfying the conditions of the above theorem, and for a target 
function / having smoothness a < 2m in Lp there is Sf~ € S'(fc,S) so that 

\\f-SfM\p = '3ih-). 

In this paper, we develop an approximation scheme delivering novel error 
estimates for a robust family of SBFs: the 'polyharmonic' kernels. This is 
the family of Green's functions of iterated and perturbed Laplace -Beltrami 
operators (see Definition 13.11 for a precise definition). Such kernels have been 
studied by Preeden and his collaborators, cf. [3j and references therein. They 
include the Green's functions for A"*, and, thus, are direct generalizations 
of the periodic "Bernoulli splines" (famously studied in [3]) and are, in 
some sense, the spherical analogues of the "surface splines" used in M'^. On 
the other hand, the SBFs obtained by directly restricting the W^'^'^ surface 
splines to S'^ are, perhaps surprisingly, often represented in this family. 

The scheme developed in this article is based on replacing the kernel 
in an integral identity by a linear combination of (few) scattered rotations 
of the kernel. This method has recently been introduced by DeVore and 
Ron in [2\ where it was used to obtain nonlinear and local results for RBF 
approximation in the boundary- free, Euclidean setting. Later, it was used 
ill [6j, [5j to provide precise approximation orders for RBF approximation 
in domains in M^. 

The layout of this article is as follows. In Section 2 we discuss some 
basics of analysis on spheres. Section 3 introduces the kernels used in this 
paper and shows that they can be expressed as a sum of surface splines. 
In Section 4 we establish a basic strategy for exchanging the kernel by a 
linear combination of its copies. Section 5 estimates the error in making 
this exchange, while Section 6 collects our main results. 

2. Background 

We denote by S'^ the unit sphere in R'^^^, and by we denote its vol- 
ume. The distance between two points, x and a, on the sphere is written 
dist(x,a) := arccos(x ■ a). The basic neighborhood is the spherical 'cap' 
C{a,p) := {x € S'^ : dist(a;,a) < p}. Throughout this article, H is assumed 
to be a finite subset of S"^, and the 'fill distance', 

h := h{E,S'^) := maxdist(Q, H), 
measures the density of H in S'^ 
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The spherical harmonic, as studied in is the basic tool of Fourier 
analysis on the sphere. For each eigenvalue, Vi := l{i + d— 1) of the Laplace 
- Beltrami operator A on S'', there corresponds an eigenspace of 'spheri- 
cal harmonics' of exact degree i, called JC^, having dimension N{d,i) := 
'''^^~^r(£+i)r[d)~^^ with orthonormal (in the sense of L2) basis {Ym/)^^^f^ ■ 
The space of spherical harmonics of degree less than or equal to L is de- 
noted Ul = Y.e<L^e- 

In this article, our focus is on zonal kernels. These are kernels on the 
sphere having the form (x, a) (p{x ■ a), where 4> : [—1, 1] — M. Such ker- 
nels, being the composition of an inner product with a univariate function, 
can be expressed in terms of an expansion in orthogonal polynomials. The 
Gegenbauer (or ultraspherical) polynomials, {Pi^^)'^q, are orthogonal on 
[— 1, 1] with respect to the weight (1 — t'^)^~^/'^ . We expand zonal functions 
on S*^ using {Pi'^'^^)'^q, with := Gegenbauer coefficients are 

1 

and the expansion is (j){x ■ a) = X^^g (^ePe^'^\^ ' Q^)- This can be expressed, 
via the addition theorem for spherical harmonics \13\ Theorem 2], as: 

00 00 N{d,e) 

e=o ^'^ £=0 m=l 

where we supplant the Gegenbauer coefficient by the Fourier coefficient 
0(£) := xf^*^^- "CLoie that the polynomials used by Miiller, [13], which 
he calls Legendre polynomials and denotes by 7^ (suppressing the depen- 
dence on d), are normalized in Loo'- they satisfy || ^P^Hl^i-i,!] =?£(!) = 1. 
The Gegenbauer polynomials used here are normalized in L2([— 1,1]; (1 — 
i2^A-i/2-j^ and are related to Miiller's Legendre polynomials by: Pj-'^''^ = 
^+2\i-i-\^p^^ Basics of Gegenbauer polynomials can be found in [151 Section 
4.7]. A key result relates the smoothness of the kernel 4> with the decay of 
its Fourier coefficients, (/>(^). 

Proposition 2.1. // E£o |0(^)K"'+^''"^ < 00 then (p G C^[-l, 1]. 
Proof. From [151 Equation (4.7.14)], observe that the derivative of a Gegen- 
bauer polynomial satisfies jiP^i {t) = 2\dP^^''+^\t). Hence, for k < i, 



-^Pl^''\t) — 2^\ci\d+2 ■ ■ ■ ^d+2k-2P£^'i^^''\t), 



while for i < k, the polynomial pj;'^'^^ is of degree at most k — 1 and is 



annihilated by ^.) Since oJdj^2^d+2 = T^^d, it follows that 

\XdU)d J at'' \uJd+2k-2 
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Utilizing a uniform bound on Gegenbauer polynomials ( ^15i Theorem 7.33.1]), 

max_i<(<i \P^^\t)\ = {^~^'^^~^), when A > we see that 

The result follows because the series 

E£o is absolutely 

convergent, and, hence, equals -^(j){t). □ 

When 1 < p < 00, the smoothness spaces we consider are the Sobolev 
(for integer smoothness) and Besov classes which we denote by VK^ (S*^) and 
-Bp oo(S'^), respectively. For p = oo,we consider C''(§'^) and the Besov classes 
-^oooo(^'^)- These can be defined on S"^ in several, equivalent, customary 
ways. The simplest way to define Sobolev spaces is to use a partition of 
unity and local changes of variables to import the definition from as in 
[101 Sect. 3]. See the reference [16] for this and other definitions. Of principal 
importance to us is the fact that Lp(S^) = W^{S'^) and A boundedly maps 
W^{§'^) to W^-'^{§,'^) (for s > 2). We postpone the discussion of Besov 
spaces until Section 6. 

3. POLYHARMONIC KERNELS AND SURFACE SpLINES 

The kernels we introduce in this section, the polyharmonic kernels, are 
fundamental solutions for certain elementary partial differential operators. 
In Section 3.1, we begin by defining the kernels in terms of the operators they 
invert. This indirect approach is taken because it is key to understanding 
the approximation scheme discussed in subsequent sections. A more direct 
expression in terms of Gegenbauer polynomials, ([3]), is also given. 

Lemma [331 provides an asymptotic expansion Gm ~ X^jt=o ^j't's+j of poly- 
harmonic kernels in terms of simpler kernels, called surface splines. This is 
developed in Section 3.3. In the course of demonstrating the asymptotic 
expansion, we make the complementary observation. Lemma 13.41 that the 
surface splines are polyharmonic kernels. This is the focus of Section 3.2. 

3.1. Polyharmonic Kernels. 

Definition 3.1 (Polyharmonic Kernels). Let m > d/2 he an integer. For 
ri, . . . , € C, the polyharmonic kernel Gm = (S'( • ; n, . . . , r^), defined on 
[—1,1), is the fundamental solution for the product of perturbed Laplace - 
Beltrami operators (A — ri) ... (A — r^). 

Our interest in polyharmonic kernels stems from certain integral iden- 
tities they satisfy. Such identities may hold for a general kernel fcjvf (not 
necessarily polyharmonic, or even zonal), 

(1) f{x)=l fi]\i{f -pf){a)kM{x-a)da + pf{x). 
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where XL^f is a differential operator of order M whose nullspace is contained 
in the finite dimensional space Ilg = X^jgg^j of spherical harmonics of 

prescribed degrees j G 3, and where pf = J2jegYlT=i\f^^j,m)Yj,m is the 
(L2) orthogonal projection onto this space. Because Hg is finite dimensional, 
< const 3) II /lip for any norm || • ||x- Hence, the identity ([T]) 
extends, by continuity, to every space T^p*^(S'^), with l<p<oo. 

Definition 3.2. // |2P holds for all f € C^'^(S'^), then kM is said to satisfy 
an integral identity of order M . 

When k2m = G{-;ri, . . . ,rm), the operator is £;2m = (A-ri) . . . (A-rm), 
and d must at least capture the indices corresponding to the eigenvalues used 
to construct ^2m- That is, 3 contains each index j G N for which there is r£ 
of the form = j{j + d—l) (there will be at most m such indices, although 
the set 8 is free to contain more). Thus every kernel G(-; ri, . . . , r^) satisfies 
an integral identity of order 2m. 

We now show that the polyharmonic kernels can be decomposed as linear 
combinations of surface splines (perhaps more accurately called "restricted 
surface splines"), which are zonal functions 



ty iog(i - 1) 



for s G N; 
s G N- i 



Roughly, these are restrictions to the sphere of a well known family of RBFs: 
the surface splines, | • |^ and | • |^log | • |, produce the fundamental solution 
of the (/3 + d)/2-fold Laplacian in R"^. The zonal kernels considered here are 
restrictions of such to the sphere, by way of the identity ^{x — a)"^ = 1 — x-a. 
Providing this decomposition is important to determining error estimates, 
because there are precise bounds for the surface splines and their derivatives, 
especially near the singularity x = a. For t > it is not difficult to see that 
there exist constants Ps,j so that 



(2) 



#(t) =/3,,^.(l-t)-^- foTj>S. 



Our investigation of polyharmonic kernels begins with observing their 
expansions in Gegenbauer polynomials. The series expansion for Gm follows 
by Fourier inversion; its Fourier coefficients are obtained by reciprocating 
the symbol of the differential operator that Gm inverts: 



Gn 



• «) = E ^ n[^(^ + d-l)- r,]-'P^'^\x ■ a). 



=1 J 



It is often useful to adopt the notation 
the Gegenbauer expansion becomes 



(3) 



Gm{x 



a] 



E 

£=1 



Hi? 



I {I, d) := £ + Arf, in which case 



Pf''\x-a). 
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3.2. Surface Splines. The series expansion for surface splines is more dif- 
ficult. It has been studied recently in ^ and [H]. These results allow a 
precise expansion of the kernel in Gegenbauer polynomials. 

Lemma 3.3. For s G N/2 satisfying m := s + d/2 € N, and i = £ + Xd, 
there is a nonzero constant Cg ( depending on s and d) such that the Fourier 
coefficient is 



m ^ 



u=l 

for i > s when d is even, and for all i when d is odd. 
Proof. The formula (j)s{x ■ a) = Yle^o Pf^'^\x ■ a) holds with 

^ l + Xd T{1 - s) 

' iOdXd r{s + i + d) 

for £ > s when s G N by [I] [(2.20)] and for all i when s G N - 1/2 by 
[1][(2.12)]. Utilizing the notation £ = i + Xd and s = s + Xd (and noting 
that s is in N — 1/2, since we assume that s + d/2 is an integer), the factor 
T{i- s)/T{s + i + d) simplifies to [{i ^ - {^f )■■■{£ ^ - s^)]~\ and the 
lemma follows. □ 

Thus, for any positive half- integer s, we have the expansion for surface 

splines: ^,(:r-a) = p{x-a) + Csj:T=olJ^, UZii^' -i'^-k)T'Pi''\x.a), 
although the extra polynomial term p G ns[— 1, 1] is only needed when d is 
even. 

Lemma 3.4. Let m = s + d/2. The kernel {x,a) ^ (f>six ■ a) satisfies an 
integral identity of order 2m with operator 

m 

= n - (i - d/W + d/2 - 1)] , 

andpf = J2i<sYlm=i\f^^i,ni)Yi,m, the projection onto Ilsm-d- 

Proof. Since the symbol of the Laplacian is ^ ^ — A^, the factor in the 
denominator of the Gegenbauer coefficient of (f)s is 

^"*^-(--^)^ = ^"'^-Al-[(.-^)^-Ai]=P-A^-(.-^)(. + ^-l). 

Thus, when d/2 is fractional, (f)s is the fundamental solution for the invertible 
differential operator whose symbol is YYjLi (^ ^ ~ ^"d) ~ (j ~ + f ~ 1 
since the eigenvalues of A are integers (the integers k{k + d — 1)). When 
d/2 is integral, the differential operator is invertible on the complement of 
the space of spherical harmonics of degree less than or equal to s. □ 



8 T. HANGELBROEK 

3.3. Surface Spline Expansion of Polyharmonic Kernels. 

Lemma 3.5. For positive integers m and d, let s = m — d/2. The polyhar- 
monic kernel {x, a) i-^ Gm{x • a; ri, . . . , rm) can he written as 

J-i 

Gm{x • ") = X] Ij 4>s+j{x ■ a) + Rj{x ■ a) 

j=0 

with Rj G C(-^+^^^)([-l,l]). 

Proof. We begin by expanding each of the Fourier coefficients of Gm- From 
dSD we observe that GZ{i) = UT=ii^^ " " ^i)"^- Factoring f-^™, we 
have, for i > max(y^|ri|, . . . , \/|rv^ , that 



j=i V / 



2m / 1 I \ ^ A /T— 2n 




The second equahty fohow by writing each factor in the product as a Neu- 
mann series (i.e., a series of the form (1 — a)~^ = Xlj^o'^"')' then by 
multiplying the m series. We do likewise for the coefficients of (ps+j (deter- 
mined in Lemma 13.31) when £ > s + J: 



r-2n 

J" 

n=l 



This allows us to choose the coefficients 70,71, .. . in succession, via 7^ := 
C^_^j{Aj — Ylk=<>^>'^j,k)- With this choice, the first J terms in the as- 
ymptotic expansion of Rj{t) are forced to vanish. The fact that each 
7j depends only on the previous coefficients 7^, k < j, is evident from 
Table [TJ The coefficients of the remainder term are determined to be 

GZ{i) = i-^^'i 1 + Ai e-^+ A2 A3 i-^ + ...) 

Mi) = Csi-^'^l 1 + Bifll-^ + B2,oi-^ + i?3/-' + • • • ) 

<gl(^) = C.+l/-2™( + 52,1^"*"' + ^3,1^-' + . . . ) 

c/)s+2{e) = Gs+2(-^H + i?3,2^-^ + . . . ) 

Table 1. Expansion of Gegenbauer coefficients 



Rj{£) = G^{l)-T.i=ll3- and 



I^J(^)I<E 



J-i 



Aj+k - E ljBj+k-j-l,j 
j=0 



£-2(m+J+fc) < ^Q^g^ (J) ^ 



00 



fc=0 

for sufficiently large £ . By Proposition 12. 1^ the lemma follows. □ 
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4. Replacing the Kernels I: Finding the Coeffigients 

We now wish to investigate a 'coefficient kernel' a : H x S'^ ^ M that wih 
allow us to effectively replace a fc(x • a) with X^^gs ^(?' '^)^(^ ' 
representation ([1]). To do so, the exchange Cfc given by: 



must be appropriately small in Loo, and it must decay away from a = x. 
The remarkable thing is that this can be achieved using only a fixed number 
of centers near to the singularity. In this section, we develop a technique 
for choosing coefficients a(^, a) that - in the following section - is shown to 
provide an appropriately small exchange. 

The two key quantities we need to resolve are the spherical harmonic 
precision (the degree of spherical harmonics reproduced by the coefficient 
kernel) and the rate of decay of the error as dist(x, a) increases. As in 
the Euclidean setting, these are related: the higher the degree of spherical 
harmonic precision, the more rapidly the exchange decays away from the 
singularity. 

Definition 4.1 (CKC). For a set of centers E C S'^ the kernel a : ExS^ ^ R 
satisfies the Coefficient Kernel Conditions (or CKC) with precision L, radius 
p and stability K if it is measurable and the following three conditions hold: 

CKC 1 (Support): a(^,a) = when dist(^,a) > p. 

CKC 2 (Precision): For S G Ul, E^gh «(^' = Sia) 

CKC 3 (Stability): max^ggd X^ges "il <K. 

Such a local reproduction property always holds for sufficiently dense 
centers. This is demonstrated in the following lemma. 

Lemma 4.2. Civen a precision L and centers H having density h < Hq 
(with ho determined by L), there exists a coefficient kernel a : H x S'^ ^ M 
satisfying the CKC with radius p = ASLp'h and stability K = 2. 

Proof. Let Hq, := H fl C(a, p), the set of centers a distance p from a. Follow- 
ing what is, by now, a fairly standard technique in scattered data approx- 
imation (originally developed for the sphere in [9j, and deftly exposited in 
|191 Ch. 3]), a coefficient kernel is shown to exist if the sampling operator 



is boundedly invertible when the domain and range are endowed with the 
Loo and ^oo topologies, respectively. To be precise, we must show that the 
norm of the inverse of the sampling operator is bounded by 2: ||-Rq^|| < 2, 
which is accomplished in Lemma 14.31 below. Bounded invertibility of Ra 
implies that the norm of the adjoint 
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is similarly bounded. By the Hahn-Banach theorem, there is a norm- 
bounded extension of the functional (Ra^) 6a € ((11^)1^^) in the space, 
(^oo("a))'- This can be viewed as an element of and, by zero exten- 

sion, it is in ^i(H). We call this sequence a{-,a) and note that its £i norm 

is bounded by || (-R~^)'|| \\6a \\ < 2. 

The measurability of the kernel a is a consequence of its piecewise conti- 
nuity, which we now demonstrate. For each i; C H, we define the open set 
Qy := ngg„C(^,/9). These can be refined to a (finite) collection of sets 

:= \ y 

that partitions S*^. For each a in Q^, the sampling operators Ra share a 
common target n|^, and the operator valued map a i— > Ra is well defined 
and Lipschitz. Indeed, \\RaP — Ra'P\\ < Cja — a'| || Vp||oo < C'lIq^ — a'| ||p||loo 
implies that \\Ra — Ra'W ^ Cilct — a'\. The inverse is similarly Lipschitz, 
because R~^—R~} = R~^ \Ra' —Ra\ R^} ■ For a fl^ the family of sequences 
a{-,a) have their support in v. To show that a ^ a(-,a) is continuous we 
simply observe that 

\\a{;a)-a{;a')\U, < \\R-^ - Ra>\\\K\\ + \\{Ra')'\\¥a' - Sa\\ 

= 0{\a-a'\). 

□ 

Lemma 4.3. Given a precision L and centers H with density h < ho, let 
p = and Hq, := H fl C(a,p) for each a € S*^. The sampling operator 

Ra is boundedly invertible on the space of spherical harmonics of degree L 
or less, and 

\\p\\L^iCia,p)) < '^\\RaP\\i^{Ec,}- 

Proof. This is accomplished by noting that spherical harmonics, when re- 
stricted to great circles, are trigonometric polynomials. From this we can 
apply the Markov inequality of Videnskii[17j, which states that for a trigono- 
metric polynomial, r of degree n 

\T'{e)\ < 2n2cot(u;/2)||r||i^(_^,^) for uJ<Tr,\e\<uj 

to control the size of a spherical harmonics having many zeros in a spherical 
cap. 

Select p G Ul and find xq such that |p(xo)| = ||p||Loo(c(a,p))- Following 
Wendland [HI p. 30], we take ^ € HnC(a, p) so that ^ is in a cone with vertex 
Xq and distance from xq less than h + < 3h (this is possible because a 
cap of radius h with center located at a distance from xq is contained 
in the cone of aperture 9). Let xq be the terminal point of the geodesic 
segment starting at xq, passing through ^ and having length p. Restricting 
p to this geodesic gives a trigonometric polynomial of degree L. Vis., there 
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is q : [— vr,7r] C, q{9) = J2\j\<2m^j^^''^ ' such that q{—p/2) = p{xo) and 
q{p/2) = p{xo) and 

^-p/2+|xo-e| 

\p{xo)-pm < / \q'{t)\dt. 



l-p/2 

By Videnskii's Markov inequahty, \q'{t)\ < cot(/9/4)||(7||^^(_p/2,p/2)) and, 
consequently, we have that 

\pixo) -piC)\ < 3h 2{Lf -\\p\\L^(c{a,p)) < l\\p\\LUCia,p))- 

Thus \\p\\^ (c{o p)) — ^11^*1=11^00 ^'^d the lemma is proved. □ 



A consequence of the CKC is that for any x and any zonal function k that 
is smooth on the interval 3x 

:= [min Qa;., max Qa;], where Qx = {x-a}U{x-S, : 
^ G (E (1 C{a, p)}, the exchange can be estimated in terms of the length of 
the interval 3x and the size of derivatives of k purely on 3x- This is the point 
of the following lemma: 

Lemma 4.4. Given a coefficient kernel satisfying the CKC with precision 
L, if k & C^^~^^\3x), then the exchange satisfies 

(4) ek(x,a) < max |x • (a - 01^+^ max Ifc^^+^^mi 

L I ^GHnC(a,p) te3^ 

Proof. Let both x and a be fixed, set x • a = G [— 1, 1] and choose the 
Taylor polynomial of degree L, qL,tai of ^ expanded about ta- Now qL,ta 
may be rewritten as a linear combination of Gegenbauer polynomials. 



£=0 



Note, furthermore, that qL,ta{x ■ a) — Yl,^{^T^)lL,tc,{x ■ i) = ^ the ad- 
dition theorem, since qL,t{x ■ () = Yl£QL~ti^)Hm^i,rn{x)Yi^rn{C) and each 
Yi^miC) is annihilated hy p = 6a — a{(,, a)6^. Consequently, ek{x,a) < 
a)| \k{x ■ ^) — qL,ta{x ■ ^)|, and the Taylor's theorem gives: 

\Hk)-QL,tAk)\<^\k-to.\'^^' max \k(^+'\u)\. 

□ 



5. Replacing the Kernels IL Estimates 

Having found coefficients suitable for replacing the kernel in a repre- 
sentation ([1]), we now obtain estimates on the exchange in an effort to 
estimate the norm of the operator '■ Lp(S'^) — > Lp(S'^), defined by 
Ek9{x) = J ek{x,a)g{a)da. The bound, \\Ek\\ := supo^ggi^ ||^fc5f||p/||c/||p 
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gives us essentially the error estimates we desire, since the pointwise error, 
£(x) := - /gd J2 a)k{x ■ ^)ii2mf{a)da\, satisfies, by ((U, 



£(x) 



< 



efc(3;,a)|iL2m/(a)|da. 



In other words, ||£||p < ||-E'fc|£'2m/| ||p < H-E'fcll ||'^2m/||p- Because of the ex- 
pansion from Lemma 13.51 we focus on obtaining the estimates for surface 
splines first, before moving to polyharmonic functions in general. 

Lemma 5.1. Let m = s + d/2. Assume a is a coefficient kernel satisfying 
the CKC with radius p, precision 2m and stability K . Then for j G N the 
exchange of the kernel (ps+j satisfies 



e<f>,+i (x, a) < const (K, m, j, d) p 



2{m+j)-d 



1 + 



dist(a;, a) 
P 



2j-d-l 



Proof. We consider three regions, for a fixed 'north pole' a € S'^: 

Qi := {x \ tt/2 < dist(x, a) < vr}, where the surface spline is smooth; 
0,2 '■= {x \ < dist(x, a) < p}, a cap of radius p near the north pole; 
^3 ■= {x \ p < dist(x,a) < it/2}: a band where high order deriva- 
tives decay. 

fii : We note that outside the spherical cap (7(Q,7r/2) the (2m + 1)*^* 
derivatives of (ps+j are bounded by /3s+j,2m+i; and we can use (jl|) to obtain 
e</,,+,(x,a) < p'^"'+^^j3s+j,2m+i- 

ft2 : In the cap nearest to a, we use the fact that (ps+j has a high order 
zero. Here we need a relationship (used later, as well) between the geodesic 
distance of two points and their inner product 



(5) l-^(dist(x,C))'<x-C<l 



(dist(x, C))^ for dist(x, C) < 7r/2. 



For even d, the proof is complicated by the log factor, so we consider this 
case only, as the odd case follows by a similar but much simpler argument. 
We proceed by writing 



ii - x-cr+^'iog|i - x-ci 



p 



l-x-C 



p- 



s+j 



log \p \ + p 



1-x-C 



p- 



s+j 



log 



1-X-C 



p- 



Since s is even when d is even, the first term is simply a spherical harmonic 
in C (by the addition theorem), of degree s + j < 2m, and is therefore 
annihilated by the functional fi = 5a — J2^eE ot(C, a)(J^. Thus, we need only 
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apply n to the second term; we obtain (by the left hand side of ([5])), 



^^^^+Aoo,a) < (1 + max 



l-x-C 



p2 



s+j 

log 



p2 



max Itl*^"' log \t\ 
0<t<9/2 



ria : The estimate ([5]) bounds the derivatives of (ps+j, but in the northern 
hemisphere, we can achieve better estimates for |x • — a)|, in the sense 
that this inner product becomes considerably smaller than p when x and a 
are close. Decompose C{a, k) \ C{a,2p) en annuli, and note that 

(6) dist{x,a) <2^p =^ |x • (a - C)| < 2'^'+^^- 

For 2''"V < dist(x,a) < 2^p we estimate maxJ^ < 1 - :^{2''-^p)'^ by © 
to obtain bounds on the (2m + 1)*'^ derivatives of (hg+i on J^- On the other 



hand, we can apply ([6]) to estimate |x • (a — ^)|. Thus, 
= const {K, m, j, d) (2'') "^'^^^^ p2s+2j 



s+j — 2771— 1 



The first inequality follows from ^ using ([2]), ([5]) and ([6]) while the second 
inequality is a consequence of the fact that d = 2m — 2s. □ 

6. Main Results 

We are now in a position to prove our main results, that polyharmonic ker- 
nels and surface splines deliver Lp approximation orders commensurate with 
the Lp smoothness of the target function, at least up to a (putative) 'satura- 
tion order': the order of the differential operator that the kernel inverts. We 
begin by giving 'high order' results, for functions of 'full' smoothness. Af- 
terwards, we give the lower orders and the corresponding smoothness spaces 
by means of real interpolation. 

Theorem 6.1. Assume the coefficient kernel a : H x S'^ — > M satisfies CKC 
with radius p, precision 2m and stability K. Assume, further, that the kernel 
k provides an integral identity ([ip of order 2m and can be decomposed as 
k = Yl'^=o^ Ij'i's+j + R, with s = m — d/2 and remainder R E C^^™") [—1, 1] . 
Then for / € H p 
approximant 



Then for f G W^"'{§'^), ifl<p<oo, or for f e C2™(S'^) when p = oo, the 



T^f{x)=pf{x) + Y,A^Hx-0, 
?e2 



with coefficients A^ = J^d^2mif ~ Pf)i<^) a('^,Q;)da, ^ G converges to f 
in Lp(S°') with error: 

11/ -7^h/||l,(s<*)< const (K,fc)p2m, 
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and with coefficients satisfying \\A\\£^(^^^ < const {K, k 

The decomposition, k = Yl'j=o^ Ij^s+j + R, means that this result holds 
for surface splines themselves, and by Lemma 13.51 it holds for polyharmonic 
kernels Gm as well. 

Proof. We begin by estimating the operator norm of E},. To do this for 
1 < p < cxD, we simply find estimates for and HE'fcHoo^oo; obtaining 

the ll-Efcllp^p norm by interpolation. By symmetry, both the Li and L^o 
operator norms are bounded by sup^g§d Jg^ |efc(a;, a)|dx. 

The decomposition of the kernel permits us to estimate this integral as 
the sum of the constituent integrals / \^<j)s+j (^j a)|da;, for j = . . . 2m — s 
and J [e/j(x, a)[dx. The latter can be estimated using Lemma 14.41 directly: 
!x\^R{x,a)\dx < 2^i^dP^™'||-R^^™'''||Loo[-i,i] • The integrals of the kernels 
^4>s+j{^^(^) ^'^^ estimated by splitting the sphere into the southern hemi- 
sphere, r^i, and northern hemisphere, Jlf . By Lemma [5. 11 e.^^^. is bounded 
uniformly over Q.i by Cp^'^^^, so /^^ \e^_.^^{x,a)\da < Cp^™+^. 

On we integrate using polar coordinates, obtaining: 

[x,a)\da < C [ p'(-+^^-'^(l + ^^^^Y~'~\a 
Jni \ P J 

( /■7r/(2p) \ 

< C7p2(m+,) / ^ i?2,-2^i2j 

To bound the coefficients, we make the estimate 

El^^l^E / |a(a,0ll^2m(/-J>/)(a)|da. 

This is less than 

E 

?6S 



□ 



a(a,0||£.2m(/ -p/)(a)|da < K \^2m{f - Pf)ia)\da 

< CK\\f\\w2m. 



The previous theorem requires the target function to have 2m derivatives 
in Lp, which is quite restrictive. To treat more general functions, we can 
first approximate a target function of lower smoothness by a nearby member, 
g, of Wp"^, and apply the theorem to g instead of /. This is an old trick 
in approximation theory, and it is a consequence of the fact that the Besov 
spaces are interpolation spaces of Sobolev spaces. We make use of the Besov 
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spaces -Bp^ooi ^ < p < oo and < a < 2m, which are the spaces of Lp 
functions with norm 

■■= s^uP (t-^ inf {ll/ - g\\p + t\\g\\w2^ : g G W^"'{S'')]) . 

When p = oo, the norm can be rewritten with (7^™" replacing Wp"^. Rather 
than paraphrase the theory here, we point the interested reader to [16| 
Chapters 1 and 7] for the pertinent theorems and definitions. 

Corollary 6.2. In the setting of the previous theorem, if f ^ Bp^CB'^) for 
1 < p < oo with < (T < 2m then dist{f , S{k, 3)) p < const p'^H/ll^a^, and 
this can be accomplished with an approximant 

with p £ Hg and with coefficients satisfying ||yl||£j(=) < Cp'^~'^"^\\f\\^a (gdy 

Proof. By real interpolation, we have, for every t > 0, that inf{|[/ — g\\p + 
t\\g\\w2^ : g G Wp2'"(S'^)} < t^-\\f\\B-^{sd)- This implies, taking t = p^"", 
that we can find gp G Wp"^{S'^) satisfying 

\\f-gp\\p < 2p'^||/b.^(s.); 

Applying the the previous theorem to gp gives ||/ — ssg\\p < 11/ — 5 Hp + 
\\g - s^g\\p < 2p-||/b.^(s.) + const p^^p'^-^^^llfW^.^^^^dy The coefficient 
estimate follows by a similar argument. □ 

By Lemma 14.21 we can apply the previous results to approximation with 
sufficiently dense centers. 

Corollary 6.3. For centers H G S'^, having fill distance h < ho, (with ho 
given by Lemma \4^ , if f £ Xp 

dist(/,5(fc,H))p< const ik)h'\\f\\x^ 

where Xp is Wp"^ when s = 2m, or Bp^^ when < s < 2m. 
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